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Effects of the quark field on the ghost propagator of the lattice Landau gauge are investigated 
by using the unquenched SU(3) configurations produced by the MILC collaboration and compared 
with quenched gauge configurations of SU(2) first copy and the parallel tempering (PT) gauge fixed 
samples and quenched SU(3) 56^ configurations. We measure the color symmetric and the color 
antisymmetric ghost propagator and the Binder cumulant of the norm and the P norm of color 
antisymmetric ghost propagators and investigate deviation from those of Gaussian distributions. 

In the first copy samples of quenched SU(2) we observe a large fluctuation in the Binder cumulant 
at the lowest momentum point. This fluctuation is reduced in the PT gauge fixed samples. The 
color anti-symmetric ghost propagator of quenched SU(3) configurations depends on the lattice size 
and is small as compared to the symmetric one in the large lattice of 56*. 

The Binder cumulants of the quenched SU(2) and the unquenched SU(3) are almost consistent 
with 3-d and 8-d Gaussian distribution, respectively. In most cases they are independent of the 
momentum, but in the case of MILC/ configuration of /3imp = 7.11, the Binder cumulant has an 
anomalous momentum dependence. This anomaly is correlated with an anomaly of mass function 
M(g) of the corresponding quark wave function of bare mass mo = 27.2MeV. 

A comparison of the SU(3) unquenched configurations and quenched configurations indicates that 
the dynamical quarks have the effect of making color antisymmetric ghost propagator closer to the 
Gaussian distribution and the Kugo-Ojima color confinement parameter c closer to 1. 
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I. INTRODUCTION 

In this presentation, we would like to show our study of 
color confinement and dynamical chiral symmetry break- 
ing via unquenched lattice Landau gauge simulation us- 
ing configuration given by the MILC collaboration in 
which Asqtad Kogut-Susskind (KS) fermions are used. 

In our study of color confinement mechanism, we mea- 
sure 1) the Kugo-Ojima confiment parameter c, 2) the 
condensate in running coupling, gluon propagator and 
quark propagator and 3) the ghost condensate parameter 
V and the Binder cumulant of the color anti-symmetric 
ghost propagator. 

We are interested in the effect of quark fields on these 
variables since the confinement and the chiral symme- 
try breaking could be related with each other. We show 
our result of measurement of mass function M{q) of 
MILCc(203 X 64 lattice) and MILC/(283 x 96 lattice) 

In the Kugo-Ojima theory 2], unphysical longitudinal 
gluons and ghosts do not appear in the physical spectrum 
by the BRST quartet mechanisme. Using arrows as the 
ERST transformations and B as the Nakanishi-Lautrup 
auxiliary field, the BRST transformation in the quark 
sector can be expressed as 

ip —ipc — > 

Inclusion of fermion gives more restiction on the degrees 
of freedom of the ghost and it may change the fluctuation 
of the ghost propagator. 



Landau gauge adopted in this theory suffers from the 
Gribov problem^, i.e. the uniqueness of the gauge field 
is not guaranteed in the simple gauge fixing procedure 
and in general there are gauge copies. In the Lagrangian 
field theory, presence of gauge copies is ignored, but in 
the lattice simulation of the infrared regionjimportance 
of the selection of a unique gauge is qualified [J, [^,0. The 
complete algorithm is however not known, despite vari- 
ous proposals (smearing Q , parallel tempering (PT) Q 
etc). 

This paper is organized as follows. In sect. II, we re- 
view theories of the color confinement in infrared QCD. 
Results of Kugo-Ojima confinement parameter in lattice 
simulation are shown in sect. III. Simulation results of the 
ghost propagator, the gluon propagator, the QCD run- 
ning coupling and the quark propagator are shown in 
sect. IV, V, VI and VII, respectively. Summary and 
discussion are given in sect. VIII. 



II. COLOR CONFINEMENT 

Let us summarize developments of the study of color 
confinement which are related to our work. In 1979 
Kugo and OjimaQ] proposed a color confinement crite- 
rion, which will be explained in the following. Gribov's'] 
and Zwanziger[4i| gave condition of the confinement on 
the infrared exponents of the gluon propagator and the 
ghost propagator. These exponents define that of the 
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running coupling as{q) in MOM scheme. A measure- 
ment of as{q) via triple gluon vertex in quenched SU(3) 
lattice simulation by the Orsay group suggested infrared 
suppression and a presence of mass-dimension 2, 
condensates [8]. 

The relation between the mass-dimension 2 conden- 
sates and the Zwanziger's horizon condition generated 
by the restriction of the gauge field in the fundamental 
modular region was pointed out by the Gent groupQ. 

Although is not BRST invariant, a mixed conden- 
sates with cc is BRST invariant [Toj . The Gent group 
suggested in the Local Composite Operator (LCO) ap- 
proach that the cc condensates would manifest itself in 
the color antisymmetric ghost propagator. 

The lattice simulation of the color antisymmetric ghost 
propagator in SU(2) lattice Landau gauge was performed 
by the group of Cucchieri[ll|. 



A. Kugo-Ojima conflnement criterion 

The Kugo-Ojima confincmnt criterion is expressed by 
the parameter c defined in the two-point function 



-iq{x-y) , 
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as ^'"''(0) = -d^^c becomes 1. We adopt the SU(3) color 
matrix A normalized as trA°'K^ = 5°''. 

The parameter c is related to the renormalization fac- 
tor of the gluon sector, the ghost sector and the quark 
sector as 



~Zz 



If the finiteness of Zi is proved, divergence of Z3 is a 
sufficient condition. If Z3 vanishes in the infrared, Zi 
should have higher order 0. li Z2 is finite zf should 
vanish. 



B. Zwanziger's horizon condition 

Zwanziger proposed realization of the unique gauge by 
restricting the configuration in the fundamental modular 
region. The condition on his two-point function given 
below is equivalent to the Kugo-Ojima theory in a naive 
continuum limit. 



^e-*«(--«) (tr U-^D^-^{-D,)k 



Here e = (E, M^"^ S{U^,^dK'')) /{(TV^ - 1)!/}, where 



S{Ux.^) is the coefficient of of the covariant derivative. 



The horizon condition reads lim G^^ {q) — e 







l.h.s. 
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0, and the 
= id-l)h 



l)c- 
— 4, and it follows 



of the condition is (^^ 

where h = c and dimension d 

d . 

that /i = — > horizon condition, and thus the horizon 
condition coincides with Kugo-Ojima criterion provided 
the covariant derivative approaches the naive continuum 
limit, i.e., e/d = 1. 



III. LATTICE DATA OF THE KUGO-OJIMA 
PARAMETER 



In our lattice simulation, we adopt two types of the 
gauge field definitions: 1) \ogU type in which the link 
variable U and the gauge field are related as Ux.f_i — e^^-'^ 
and 2) [/—linear type in which A^,^ is a traceless part of 
the difference of Uxfi and C/J ^ multiplied by one half. 

The optimizing function corresponding to the two def- 
initions are 1) The P norm of A^ and 2) tr{2 - {U^x,,^ + 

U^x,fj.^)), respectively. 

Under infinitesimal gauge transformation g^^Sg = e, 
its variation reads for either defintion as 

AFuig) = -~2{dAa\e) + (e\ - dD{US)\e) + • • • , 

where the covariant derivativative D^^ {U) for two options 
reads commonly as 



X,fJ./ 



= s{Ux,^)d^^ + [Ax,^,4'] 

where d^cj) = (j){x + ^i) — 4>{x), and = 0(3: -I- -|- (j){x) 
The function S{Ux.^) in the \ogU definition is 

q/Tj \ _ 

^11^ a:/.; - tanh(adiA,,^/2) ' 

Stationality of F[j{g) means Landau gauge, the local 
minimum means the Gribov region and the global mini- 
mum means the fundamental modular regionQ. 

In the quenched Landau gauge QCD simulation, the 
Kugo-Ojima parameter saturated at about 0.8, as shown 
in TABLE [D The parameter e/d is closer to 1 in the 
log U definition of the gauge field. 

We measured the corresponding values of the un- 
quenched configurations of MILCc and MILC/[14]. The 
lattice specifications are shown in TABLE [Til 

As shown in TABLE HITl the Kugo-Ojima parameter of 
MILC configurations are consistent with 1. The qualita- 
tive difference from quenched simulations would be due 
to the difference of the ghost propagator caused by the 
quark field. 



IV. THE GHOST PROPAGATOR 



= G 



The ghost propagator is defined as the fourier trans- 
form of the matrix element of the inverse Faddeev-Popov 
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TABLE I: Kugo-Ojima parameter and Zwanziger parameter 
in L*'— Imear(left) and log (7(right). /3 = 6.0 and 6.4. 



p 


L 


ci ei/d h\ 


C2 62 /d h2 


6.0 
6.0 
6.0 


16 
24 
32 


0.576(79) 0.860(1) -0.28 
0.695(63) 0.861(1) -0.17 
0.706(39) 0.862(1) -0.15 


0.628(94) 0.943(1) -0.32 
0.774(76) 0.944(1) -0.17 
0.777(46) 0.944(1) -0.16 


6.4 
6.4 
6.4 


32 
48 
56 


0.650(39) 0.883(1) -0.23 
0.739(65) 0.884(1) -0.15(7) 
0.758(52) 0.884(1) -0.13(5) 


0.700(42) 0.953(1) -0.25 
0.793(61) 0.954(1) -0.16 
0.827(27) 0.954(1) -0.12 


6.45 


56 




0.814(89) 0.954(1) -0.14 



TABLE II: Pimp, the bare quark mass mo, the inverse lattice 
spacing 1/a, lattice size and lattice length(fm). Suffices c 
and / of MILC correspond to coarse lattice and fine lattice. 

Pimp ~ 5/3 X [3. 







amoiud) amo^s) 


l/a(GeV) Ls Lt 


aLs(fm) 


MILCc 


6.83 


0.040 0.050 


1.64 


20 64 


2.41 




6.76 


0.007 0.050 


1.64 


20 64 


2.41 


MILC/ 


7.11 


0.0124 0.031 


2.19 


28 96 


2.52 




7.09 


0.0062 0.031 


2.19 


28 96 


2.52 



operator. 



FT[D-^{x,y)] = FT{tr{k^{{M[U])-').,yA')), 

where M = —d^D^. The ghost dressing function G{q^) 
is defined as q^Dc[q^). 

We calculate the overlap to get the color diagonal ghost 
propagator 



DG{q) 



1 1 



xS^\{A'' cosq • x|/,^(x)) + (A'^ sinq • x|/,^(x))) 
and color anti-symmetric ghost propagator 
1 1 
JfV 

xr^^((A''cosq.x|/,^(x)) - (A"sinq.x|/,^(x))) 



where = 2 for SU(2) and 6 for SU(3). Here /c''(x) and 
//(x) are the solution of M fix) = p''{x) with p''{x) ^ 
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A" cos q • X and A" sin q • x, respectively. 
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FIG. 1: The ghost dressing function of MILC/ Pimp ~ 
7.09(stars) Dashed line is the 4-Ioop Nf = 2 pQCD result 
(Ag = 3.01, y = 0.0246100) 

The ghost dressing function of MILC/ Pimp — 7.09 
and that of the quenched /3 = 6.45 56*, the color sym- 
metric and anti-symmetric ghost propagator of MILCc 
are shown in fl§\. FIG [T] is a fit of the ghost dressing 
function of MILC/ by the perturbative QCD(pQCD) in 

MOM scheme 111. 



A. The ghost condensate 



In the LCO approach, a preference of a specific di- 
rection of (/"''^c^c'^) in the color space is regarded as a 
signal of the ghost condensate. We define a parameter 
V oc \{f'^'"^c''c'^)\ and following [ll[ measure its value as 
follows. We take into account the finite size effect of the 
lattice through a parameter r, defined as 



cos[Trq/L) q^ 



where L is the lattice size and q — 0,1,---L. In our 
asymmetric (L^ < Lt) lattice and the momentum q near 
diagonal in the 4-d space, we consider the sine momentum 



TABLE III: The Kugo-Ojima parameter for the polarization 
along the spacial directions Cx and that along the time direc- 
tion ct and the average c, trace divided by the dimension e/d, 
horizon function deviation h of the unquenched KS fermion 
(MILCc, MILC/). The log [/ definition of the gauge field is 
adopted. 









Ct 


c 


e/d 


h 


MILCc 


6.76 
6.83 


1.04(11) 
0.99(14) 


0.74(3) 
0.75(3) 


0.97(16) 
0.93(16) 


0.9325(1) 
0.9339(1) 


0.03(16) 
-0.00(16) 


MILC/ 


7.09 
7.11 


1.06(13) 
1.05(13) 


0.76(3) 
0.76(3) 


0.99(17) 
0.98(17) 


0.9409(1) 
0.9412(1) 


0.04(17) 
0.04(17) 



(2sm— ) 

Ijnr 



(2sinfif 



and choos e i n the numerator as the square root of the 
volume y^L^Lt 

The absolute value of the color antisymmetric ghost 
propagator is parametrized as 



1 



1 ^ 



3 



TABLE IV: The fitted parameters r,z and v of the color antisymmetric ghost propagator \4>[q)\ of MILCc and MILC/. The 
parameters f,z and v are those of the color antisymmetric ghost propagator squared 4>{qf ■ Two values of U of MILC/ 
correspond to the average below q = IGeV and the average above IGeV, respectively. 





mo(MeV) 


r z V 


r z V 


U 


6.76 
6.83 


11.5/82.2 
65.7/82.2 


37.5 3.90 0.012 
38.7 3.85 0.007 


33.5 7.6 0.045 
33.5 7.6 0.048 


0.53(5) 
0.57(4) 


7.09 
7.11 


13.6/68.0 
27.2/68.0 


134 3.83 0.026 
112 3.81 0.028 


251 7.35 0.044 
164 7.34 0.002 


0.57(4) /0.56(1) 
0.58(2)/0.52(l) 



o -2 

Cn ^ 
O ~J 

-4 



-0.7^0.^0.25 0.25 0.5 0.75 

Log^oq (GeV) 



FIG. 2: logio \(f){q)\ as the function of logio(g(GeV)) of MILC/ 
and the fit using r = 134 and v = 0.026GeV^. 



We tried to fit \(f>{q)\ and of MILC/ f3^mp = 7.09 

and 7.11 and found a relatively good parameter fit of 
|0(q)| with use of r = 134 and v — 0.026GeV^, as shown 
in FIG El 

Extraction of v from lattice data is not so easy due to 
the finite size effect of the asymmetric lattice. Except 
MILC/ Pimp = 7.09 which had difficulty in the fitted 
value from 0(q)^ is larger than that from |0((?)|. 



B. Binder cumulant 

Cucchieri et alflj measured the Binder cumulant of 
the color antisymmetric ghost propagator 



Uiq) 



1 



and found that it is about 0.44, which is compatible with 

the expectation value of the d-dimensional Gaussian dis- 

. (04) d + 2 . , , ^ 
tribution, — = with a = 3. 

The Binder cumulant U{q) of the color anti-symmetric 
ghost propagator (/)^(q) of MILC/ (3 = 7.09 mo — 
13.6MeV/68MeV and MILCc mo = 11.5MeV/82MeV are 
stable but the data of Pimp — 7.11 is rather noisy. In the 
region q > IGeV, the Binder cumulant U{q) is 0.52(1) 
but in the region q < 1 GeV it is 0.58(3), which is con- 
sistent with that of Gaussian distribution U{q) ~ 0.583. 
The anomalous momentum dependence occurs also in the 



mass function of the quark propagator, which is discussed 
in the sect. VII. 

We compared the Binder cumulant of SU(2) color anti- 
symmetric ghost propagator of samples which are gauge 
fixed by the PT method [| and the first copy. The first 
copy shows large fluctuation at the lowest momentum 
point. The Binder cumulant of quenched SU(3) seems to 
depend on the size of the lattice. In the 56'' lattice sim- 
ulation, color anti-symmetric ghost propagator is very 
weak and the Binder cumulant is noisy. 

In [l^ we showed the Binder cumulant of MILCc 
(21 samples). We increased the statistics of MILCc to 
Pimp = 6.76 and 6.83, 20 samples each and extended 
the measurement to MILC/ Pimp = 7.09 and 7.11. In 
MILCc P — 6.76, we found an exceptional sample which 
we encountered also in the quenched SU(3) 56'* lattice^. 
The ghost propagator of this sample has large exponent 
ac which is defined at about q ^ 0.4GeV. In the case of 
quenched simulation such copy had larger || than the 
average and the PT gauge fixing excludes such samples. 



V. THE GLUON PROPAGATOR 

The gluon propagator is defined as the two-point func- 
tion 



(5,. ~ ^)DA{q')S^'. 



The gluon dressing function is defined as Z{q^) = 
q'DAiq^) 

The gluon propagator and the dressing function of 
MILC/ Pimp = 7.09 are calculated in [ij]. The data 
suggest that the gluon propagator is infrared finite. The 
infrared exponent ac and it is close to —2aG- The expo- 
nent K of the Dyson-Schwinger equation(DSE) approach 
is defined at the lower momentum point, inaccessible in 
the present lattice simulation. We observe the numer- 
ical value of ac is about k/2, and the tendency that 
qd decreases but > — 1 i.e. the gluon propagator is 
not infrared vanishing. The combination of ao + ^ac 
is supposed to approach in the DSE approach, but on 
the lattice it is about -0.1 in the MILCc and -0.2 in the 
MILC/. 
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VI. THE QCD RUNNING COUPLING 

In the DSE approach, the ghost-gluon couphng in the 

M OM scheme is calculated by the gluon dressing func- 
tion Z3 and the ghost dressing function Z3 and the vertex 
renormalization factor Zi as 



g{q)^Z-'Z^^\^^^q')Z,{^,',q')g{^^) 



Measuring the running coupling from the ghost-gluon 
coupling is more stable than that of the triple gluon ver- 
tex. Our lattice simulation (T^ of the gluon propagator 
and the ghost propagator of MILCc yield the running 
coupling shown in FIG (21 There are deviation from the 
pQCD (dash-dotted line) and the DSE approach with 
K = 0.5 (long dashed line). As was done by the Orsay 
group 0, we consider a correction including the con- 
densates and obtained (A^) ^ a few GeV^. 




-0.4-0.2 0.20.40.60.8 1 

Log_10 [q (GeV) ] 

FIG. 3: The running coupling as{q) as a function of 
logiog(GeV) of MILCc (a = 0.12fni) = 6. 76 (triangles) 

and 6.83(diamonds), (50 sarnies each). 

The running coupling in the infrared can be estimated 
from the quark-gluon coupling. 

giq) ^ Zr'zl^'{^\q')Z,i^^q')g{p). 

where Z2 is the quark dressing function and Z'^i is 
the vertex renormalization factor. An evaluation of 
Z2 (/i^ , 9^ ) is given in the next section. 



VII. THE QUARK PROPAGATOR 

We extended the measurement of the quark propagator 
using Asqtad action of MILCc [H to MILC/. In the 
case of MILCc, we compared the Asqtad action and the 
Staple+Naik action. 

Due to long computation time for the convergence of 
the conjugate gradient method, the number of samples is 
of the order of 10 for each (3im,p and the bare quark mass 
mo. 

The quark propagator is defined as a statistical average 
over Landau gauge fixed samples 



Sapip) = ( {Xp,a\ 



1 



'ip{U)+mQ 



In this expression, the inversion, -j^r — , is performed 

via conjugate gradient method after preconditioning, and 
we obtain 

The mass function M{q) reflects dynamical chiral sym- 
metry breaking. In high momentum region, it is 
parametrized as 



M{q) 



m(M^)[log(^VA|^^)]'^" 
[log(gVAlcD)]''" ' 

where dM = 12/(33 — 2Nf) and m(/i^) is the running 
mass. 

In the infrared region, we adopt the monopole fit 



M{q) 



gA3 
g2 + A2 



Too 



The momentum dependence of M{q) and Z2{q) of 
Too = 13.6MeV in the infrared region of Asqtad action is 
smoother than that of the Staple-|-Naik action. It could 
be attributed to the effect of the tadpole renormalization. 
The parameters c and A in our fit of the mass function 
are given in TABLE |Vl 

We showed the quark wave function renormalization 
Z^q^) = 51 (m') 7^2(9') ofMILC/ P,mp = 7.11 using 
the staple-t-Naik action in [1^ , where Z2 is the bare lat- 
tice data and gi{q^) is the coefficient of 7^ of the vector 
current vertex that compensates the artefact in Z2. 

We adopt (A^) as a fitting parameter and calculate [1] 



Z^{q') 



2^ ^ .gi(M^) 



C2 



where a{q) are data calculated in the MOM scheme us- 
ing the same MILC/ gauge configuration ^6|]. 

Here Nf is chosen to be 2 but the data does not change 
much for 3. We choose Aqcd = 0.691GeV and (V'V')p = 

Since gi{q^) in the infrared is expected to be given by 
the running coupling, no suppression of the quark wave- 
function renormalization suggests that the infrared sup- 
pression of the running coupling obtained by the ghost- 
gluon coupling could be an artefact. 

In [1^ the ^2(9) is normalized to 1 at q = 3GeV. In our 
simulation without this kind of renormalization, Z2{q) at 
q = 3GeV is close to 1 and the results are consistent. 
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TABLE V: The parameters c and A of the Staple+Naik action(left) and the Asqtad action. 



Pimp 


mo(MeV) 


c 


A(GeV) 


cA(GeV) 


c 


A(GeV) 


cA(GeV) 


6.76 


11.5 


0.44(1) 


0.87(2) 


0.38 


0.45(1) 


0.91(2) 


0.41 




82.2 


0.30(1) 


1.45(2) 


0.43 


0.33(1) 


1.36(1) 


0.46 


6.83 


65.7 


0.33(1) 


1.28(2) 


0.42 


0.35(1) 


1.25(1) 


0.44 




82.2 


0.30(1) 


1.45(2) 


0.43 


0.33(1) 


1.34(1) 


0.45 


7.09 


13.6 


0.45(1) 


0.82(2) 


0.37 


0.50(2) 


0.79(2) 


0.39 




68.0 


0.30(1) 


1.27(4) 


0.38 


0.35(1) 


1.19(1) 


0.41 


7.11 


27.2 


0.43(1) 


0.89(2) 


0.38 


0.20(2) 


1.04(3) 


0.21 




68.0 


0.32(1) 


1.23(2) 


0.40 


0.36(1) 


1.15(1) 


0.42 



1.4 
1.2 

1 ■ 

5 . 8 
N . 6 
0.4 
0.2 



♦ 



2 3 
q[GeV] 



. 5| 



1 H 



12 3 
q[GeV] 



FIG. 4: The wave function renormalization factor ^2(9) of 
MILC/(/3i™p = 7.11 blue diamonds) and MILC/(ft™p = 7.09 
red stars). 



FIG. 5: The mass function M{q) of MILC/(/3i,„p = 7.11 
mo = 27.2MeV, blue diamonds) and MILC/(/3™p = 7.09 
mo — 13.6MeV, red stars). 



imp 



7.09 are about 20% 



Our mass function M{a) of (3,, 

larger than those of but if M{q) is renormalized to 
the theoretical value of at g = 3GeV we reproduce 
their data. Our mass function is consistent with DSE 
analysis 

In the case of /3 = 7.11 mo = 27.2MeV, the mass func- 
tion of Pimp = 7.11 Too = 68MeV is consistent with p^ . 
but that of Too = 27.2MeV is suppressed in the infrared 
and differ from [Ts]. The momentum dependence is cor- 
related with the momentum dependence of the Binder 
cumulant of color anti-symmetric ghost propagator. The 
disagreements with [l8!| may be due to Gribov copies. 

In the case of [3imp — 7.09, the chiral symmetry break- 
ing mass measured by using the Staple-|-Naik action and 
measured by using the Asqtad action (FIG [5]) are consis- 
tent within about 5%, and there is no anomaly. Whether 
the Asqtad action has a non-QCD like behavior in a cer- 
tain parameter region is under investigation. 



VIII. SUMMARY AND DISCUSSION 

In this presentation, we showed our results of un- 
quenched lattice Landau gauge QCD simulation using 
the KS fermion and infrared features of the ghost prop- 
agator, the gluon propagator and the quark propagator. 
We studied the Kugo-Ojima confinement criterion and 



. 5 
0.48 
0.46 
> 0.44 

" . 42 

< 

0.4 
.38 
.36 



40 60 
mO [MeV] 



80 



100 



FIG. 6: The chiral symmetry breaking mass cA as a function 
of bare mass and its chiral limit. Dotted line is the extrapo- 



lation of MILC/ /3i„ 



7.09 and the dash-dotted line is that 



of MILCc, Asqtad action. 



the Zwanziger's horizon condition. 

The infrared exponent of the ghost propagator ctp of 
the lattice is about half of the value k of the DSE[19I.[2oI|. 
The gluon propagator is infrared finite. The Kugo-Ojima 
parameter c was saturated at about 0.8 in the quenched 
SU(3) but became consistent with 1 in the unqucnchcd 
SU(3). 

In /3 = 6.4, 56* quenched SU(3) configurations, we 
found a copy whose = 0.272 v.s. aG(average)=0.223, 
and whose gluon propagator has an axis along which 
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the reflection positivity is manifestly violated. In the 
Pimp = 6.76, 20'^ X 64 MILCc configurations, we find also 
an exceptional sample. The exceptional samples of un- 
quenched configuration causes fluctuation in the Binder 
cumulant of color antisymmetric ghost propagator in the 
infrared, similar to the SU(2) first copy. 

In the measurement of the running coupling as{q'^) 
by the ghost-gluon coupling, we observed infrared sup- 
pression. The quark wavefunction renormalization Z^{q) 
in the infrared suggests that the infrared suppression is 
a lattice artefact and the running coupling in the con- 
tinuum limit freezes to a finite value. In the analysis, 
the artefact inherent to the compactness of the manifold 
should also be taken into account (2^ . 

In the ghost sector, the Kugo-Ojima confincnmcnt cri- 
Z 1 

terion — — = is satisfied. In the gluon sector, 
Z3 00 

the continuum limit of Z3 in the infrared is unknown, 
but the finite lattice data suggests that it is finite and 
non-zero. The infrared limit of Zi is not known, but the 
behavior of the ghost propagator suggests that it could 
be indefinite (fluctuating). 

A natural question, 'How do the fixed points of the 
QCD renormalization group flow behave in the infrared?' 
may arise. In 1971, Wilson argued that the fixed point of 
the renormalization group flow of QCD could have limit 
cycle structure rather than a fixed point ^23i]. The devi- 
ation of the ghost dressing function from 4-loop pQCD 
result in the infrared region shown in FIG[l]suggests that 
the flow of the ghost wavefunction renormalization is not 
a trivial one and it is correlated with the presence of 
dynamical quarks. The theorem on no renormalization 



effect on the ghost propagator [2^ does not guarantee the 
validity of the simple extension of the ghost wave func- 
tion renormalization from ultraviolet to infrared. 

An extension of the analysis to flnite temperature is 
straightforward. A preliminary simulation using MILC 
finite temperature configurations [2^ shows that the 
Binder cumulant at finite temperature becomes smaller 
than that of Gaussian distribution as temperature rises. 
The effect of quenching randomness [2^ played by quarks 
in the zero temperature unquenched configurations be- 
comes weak above the critical temperature Tc- To clarify 
the nature of the infrared fixed points, systematic studies 
of finite size effects and the Gribov copy effects (6l.[27l| are 
necessary. 

We observed that differences of the color antisymmet- 
ric ghost propagator in the quenched/unquenched config- 
uratiions and zero/finite temparature make differences in 
Kugo-Ojima parameter. The results imply that the quark 
field and the ghost field are correlated as the BRST quar- 
tet mechanism suggests. 
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